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When a magnetic field is applied along a direction deviated from the quantization Z-axis, the 
conservation of total magnetization holds no more. In this case the inclined field can cause a strong 
spin-evolution via the linear Zeeman term even the field is as weak as a percentage of mG. An 
approach beyond the mean field theory is proposed to study the evolution of small 87 Rb condensates 
under the weak inclined fields. The time-dependent populations of spin-components are given in 
analytical forms. The evolution is found to be highly sensitive to the magnitude and direction of 
the field. 
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II. ONE-BODY SYSTEM 
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Since the experimental realization of spinor Bose- 
Einstein condensates [l|, 0, IE EL IE IE] , the spin-evolution 
of the condensates has become a hot topic due to its 
academic interest and potential application. The manip- 
ulation of the evolution is a central problem^, H, H, EE, 
El HI El El EE EE]- Usually, one applies an exter- 
nal magnetic field lying along the direction of the axis of 
quantization (Z-axis)[g, El EE HE]- In this way the to- 
tal magnetization of the condensate is conserved, hence 
the linear Zeeman term of the field plays no role. How- 
ever, if the direction of the field deviates from the Z- 
axis, the conservation of the total magnetization will not 
exist. Accordingly, linear Zeeman term will affect the 
evolution. Since the linear term is much stronger than 
the quadratic term, a very weak inclined magnetic field 
might cause remarkable effect. The dynamic response of 
the condensate to a transversal magnetic field has been 
studied by Yi and Pu[H|22|]. Strong oscillation has been 
found in the evolution of populations of spin-components. 
Spin squeezing and macroscopic entanglement have been 
found in the studies of the ground state structures. Since 
the direction of the field might affect the dynamic phe- 
nomena of condensates sensitively, this topic deserves to 
be further studied. 

The aim of this paper is to study the effect of an in- 
clined magnetic field (lying along an arbitrary Z'-axis) 
on the spin-evolution of the condensates of 87 Rb atoms. 
Recently, more attention was paid to the study of small 
condensates with lower densities because the dissipative 
process can be reduced, and therefore richer phenomena 
might be observed[24|]. For these small systems the valid- 
ity of the mean field theory might be limited. Therefore, 
in this paper, an approach beyond the mean field theory 
is proposed. It turns out that analytical solutions can be 
obtained as follows. 



In order to understand better the effect of the inclined 
B on many-body systems, we study firstly a simple ex- 
ample, namely, the evolution of a single spin- 1 87 Rb atom 
under the field. The initial spin-state of the atom Xn 1S 
defined in a Z-frame, while B is lying along another axis 
Z' . The angle between Z and Z' is 8. The Hamiltonian 
H = —pSz 1 , where p = jB, and 7 = gpfig being the gy- 
romagnetic ratio. The quadratic Zeeman term is much 
weaker than the linear term (e.g., if B = 1 mG, the for- 
mer is 10~ 7 times weaker). Therefore, when Z' and Z do 
not overlap, the former can be neglected. 

Due to the inclined magnetic field the spin-state at 
time t becomes 
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A.\J 



(1) 



where r = pt/h, x'u is defined in the Z'-frame, and 
6) is an element of the well-known rotation matrix. 
Since the observation is made in the Z-frame, Eq. ([T]) is 
rewritten as 

= ^M^XA (2) 



where 



M^(t) = Y J dl,{-ey-d\M (3) 



P\[t) is the probability 



It is obvious that \M^\(t)\ 2 
that an atom in \x initially would be in A at t due to 
the inclined B. The evolution appearing as a variation 
of Px(t) is strictly periodic with the period r = 2tt or 



t p . For an example, when B = 1 mG, t p 



A 

t = h/p 
1.41 ms. 

It is obvious that {Mn\} is a unitary matrix and is 
symmetric. Due to the symmetry of ci* , we have 



*The corresponding author 
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In particular, 

1 + cos 2 I 
Mil = 
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+ i cos 9 sin(r) (5) 
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The other elements of {M^\} can be obtained via 
Eq. Q. Thus the evolution is completely clear. 

Incidentally, based on the mean field theory, the ef- 
fect of an inclined magnetic field has been studied by 
Pu, et al[2l|. When the atom-atom interaction and the 
quadratic Zeeman term have been neglected, they have 
derived a set of dynamic equations for the time evolution 
of the field amplitudes a\(t). When the condensate is ini- 
tially fully polarized (ai(0) = 1 and ao(0) = o_i(0) = 0), 
the set a\(t) has an analytical solution as shown by 
Eq. (0 of their paper. It turns out a\(t) = Mi.\- Thus, 
for this case, the mean field theory for many-body sys- 
tems and the above simple consideration for a single-body 
system lead to the same result. 



is seriously affected by the number of atoms initially in 
fi = 0. From Eq. (JTUJ), we know that whether N^/N is 
larger or smaller than 1/3 is crucial to the evolution. If 
= N/3, P ' would remain constant (without evo- 
lution). Otherwise, Pq' 6 will oscillate around a back- 
ground, and the amplitude would become the largest if 
Nq = N. On the other hand, if Z' and Z overlaps, we 
have X = and P ' remains also constant as expected. 

The time-dependent magnetization N(P(' e — PL\) = 
NPmag can be obtained from Eq. @, we have 

P^ g (t) = M'[cos 2 9 + sin 2 9cos(t)]/N (12) 

It implies that the time-dependent magnetization de- 
pends on the initial magnetization M 1 but is not affected 
by Nq/N. In particular, if the system is zero-polarized 
initially, Pma g (t) remains zero (without evolution). On 
the other hand, when Z' and Z overlaps, Pmag(t) remains 
to be a constant M 1 /N as expected. 

Incidentally, when 9 = ir/4, the evolution of P^' 6 (t) has 
been calculated numerically in [2l| and plotted in Fig. 3 of 
their paper. We found that the difference between their 
numerical results and those from Eq.([9|) is very small. 
The small difference implies that, in the early stage of 
evolution, the strong oscillation shown in their figure is 
essentially caused by the inclined field and is less affected 
by the interaction. 



III. MANY-BODY SYSTEMS WITHOUT 
ATOM-ATOM INTERACTION 

The initial state is assumed to be a Fock-state defined 
in the Z-frame as \I) = {N^N^Nl^, where is the 
number of atoms in \^ initially. When the atom-atom 
interaction is neglected, it is straight forward to obtain 
the probability of an atom in A as 

pi'\t)=Y,\ M ^ N 'i N ( 9 ) 

As before, the observation is made in the Z-frame. A 
notable point is that Pb does not depend on N but the 
ratio N v /N, and it has the same period h/p as the single 
atom has. Furthermore, it is invariant under a reflection 
against the X-Y plane, P^' e = P?'*~ 6 . In particular, for 
A = 

P^ =X(l-3§) + § (10) 

where 

' 2 

I = ^{ 1 - sin 2 9 cos 2 (r) + cos 2 [1 - 2 cos(r)]} (11) 

Eq. (JTOJ) implies that the probability is not at all af- 
fected by the initial magnetization N[ — Ni ± = M 1 , but 



IV. MANY-BODY SYSTEMS WITH 
ATOM-ATOM INTERACTION 

For realistic condensates of 87 Rb atoms as an example, 
the interaction 

v l3 =5{v l -v. ] )Y,9sP: j (13) 

where the strength g s — ATrh 2 a s /M, M is the mass of 
atom. P?j is the projection operator of the s-spin-channel 
(s = or 2, which is the total spin of the two atoms i 
and j). a = 101. Sag, a 2 = 100.4a_B from 23]. It is 
assumed that the number density of the condensate and 
the temperature are sufficiently low so that the single- 
spatial-mode approximation can be adopted[H. Under 
this approximation, when an irrelevant constant and the 
quadratic Zeeman term have been dropped, the Hamil- 
tonian reads [l(| 

H = GS 2 -pJ2£z>i (W) 

i 

where S is the operator of the total spin of the many- 
body system, G=\(g 2 - go) J dr\(j)(r)\ 4 , <j>(r) is the spa- 
tial normalized wave function of an atom (all atoms are 
assumed to condense into this state) . Since the details of 
4>(r) affects only the strength G, it is not essential to our 
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qualitative results. Therefore, it is simply evaluated via 
the Thomas-Fermi approximation. 

Let us introduce the total spin-state i?s,m of the whole 
system with conserved total spin S and its Z-component 
M. The overlap of this state and the initial state 
(i?S,m|-Ni,No>-W-i) = 5 MtM iD%j? N i is essential to the 
following calculation. These coefficients have been given 
explicitly by Wu (Eqs. (36) and (37) of HI). They could 
be derived also by a set of recursion formulae [2f|. On the 
other hand, we define further the total spin states i?g M , 
relative to the Z'-frame, $' s M , is related to $s,m via a 
rotation. Thereby the initial state defined in the Z-frame 
can be expanded by z9g M , as 



S,M' 



S,M' U S,M 



(15) 



where 



d n,s d s 



A- 



Since the set M , are the eigen-states of the Hamil- 



M 1 ,NJ[ u M> ,M' 
'S,M' 

tonian with the eigen-energy GS{S + l)—pM , the time- 
dependent solution of the system reads 

*(*) = E Cs^- iiGS{S+1)/p - M ' )T #kw (16) 

S,M' 

By using the fractional parentage coefficients given in 
[27l |28| , we can extract the spin-state of a single particle 



(say, particle 1) from $' s M , as 



-B(N,S,M\(i)tif_-%,_ 



(17) 



where 



_ , (N-S)(S+l) n/2 „sM' 
N(2S + 1) 



A(N,s,M',fj,) = r\ T ~^: ~' Y rz CI +1 , M ,_^ 



(18) 



sW^'d = [ %g+t) 1) ]1/a ^^ (19) 

where the Clebsch-Gordan coefficients have been intro- 
duced. Note that in Eqs. (JTSJ) and HU) N — S must be 
even, otherwise the state "d' s M , does not exist. 

Since the observation is made in the Z-frame, the single 
particle state Xu(l) m Eq. (|T7j) is further rewritten as 
= X)a u(^)XA) where xa is defined in the Z-frame. 

With these transformations, eventually the probability 
of a particle in A can be extracted from \& (t), and we have 

^*) = E<^)4/^) E ^"-M',M'-, 

S,M" ,M' 

{Cl%C I s %[A{NSM' l ^)A(NSM'^) 
+B{NSM"n / )B(NSM , t i)} cos((M" - M')r) 
+ 2C 's+2,M» C s;L" B ( iV ' 5 + 2 ' M"(i')A(NSM'fj,) 
cos([4G(5 + 3/2)/p- (Af" - Af')]r) } (20) 



0.5 



0.4- 



P 



0.1 



0.0* 
0.1 



p le 0.0 



-0.1 



■ /T 


\ / / \ \ 


■ / '■ 


1 1 1 

\ //w 


/ 

- f 


M\ 

\ K\ 


A 

// / 


>H X 


| / 

1/ 


\ \ 
V 


/ 

1 // 


/ V 


7/ 

// 1 


\ \ 


\\ / / 


I 4 / (a) 


— i — i — i — i — i — i — i — i — i — i — 


1 






2 // 

4. 5,-'' 
— I 1 1 h 


\ 

1 1 


— 

(b) - 



50 



100 150 

t(ms) 



200 



250 



FIG. 1: (Color online) P'{° (t) (a) and P^ e ag (t) = P/' 9 - Pj-f 
(b) against t calculated from Eq. (|20|l with uo = 300 x 2tt, 
N = 100 and B = O.OlmG. The initial state has M 7 /iV = 0.1 
and N/, = 0. For the curves "1" to "5", 6 = 0, tt/6, it/3, tt/2, 
and tt/2, respectively. P^ ,7T ~ e (t) = P^' e (t) holds always. The 
curve "1" is identical to the one with B — 0. "5" is for the 
case with the interaction ignored (G = 0). 



This is a generalized version of Eq. ([9]) with the realistic 
interaction taken into account. It is also a generalization 
of Eq. (3) of [29( , taking the effect of the inclined mag- 
netic field into account. When 8 = 0, Eq. (|2"D|) is identical 
to Eq. (3) of [H. 

Eq. ([20| provides an analytical description of the evo- 
lution and all the coefficients involved have analytical 
forms. Based on Eq. (|20[) . numerical results are shown 
by the following figures as examples to demonstrate the 
feature of evolution. The condensate of 87 Rb atoms is 
assumed to be trapped by a harmonic potential with 
lu = 300 x 2tt. A very weak field with B = 0.01 mG 
is chosen. This is enough to show the effect of the linear 
Zeeman term. Although Eq. ([20]) holds for arbitrary N, 
numerical results are limited by the ability of computer. 
As the first example, N = 100 is chosen. 

In Figs. [T] and [2 8 is given at a set of values. Related 
to the figures the following points are noted. 

(i) Comparing curves "2" to "4" (in solid line) with 
"1" (in dash line), the strong and swift oscillation caused 
by the inclined weak magnetic field is impressive. The 
evolution is highly sensitive to 8. 

(ii) When the interaction is ignored, the evolution is 
described by the curve "5" (in dotted line) which arises 
purely from the inclined field. Comparing "4" and "5" 
(both have 8 = tt/2), we know that the effect of inter- 
action is weak in the early stage because the two curves 
almost overlap. However, if t is larger, the influence of in- 
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FIG. 2: (Color online) All the same as Fig.[TJi but with Nq = 
0.9. 

teraction would become more and more explicit, and the 
deviation between "4" and "5" would be serious (refer to 
Fig.©. 

(iii) The evolutions shown by the curves (except "1") 
are nearly periodic with the period t p = h/p = 141 ms. If 
the interaction is removed or if G/p is an integer, it would 
be exactly periodic as shown by Eq. (f^Uj) . However, since 
t p is reversely proportional to B, the cycles of oscillation 
would become very dense if B is large. For instance, 
when B > 1 mG, the curves of P x ' (t) with 8 ^ look 
like a band. The width of the band is equal to a double 
of the amplitude of oscillation, and therefore depends 
on 8. However, if 8 is not small (say, 8 > 7r/15), the 
width would be too broad and the band is difficult to 
be determined. Thus, when a stronger field is used in 
experiments, only the cases with small 8 are meaningful. 

(iv) When the field is not inclined, magnetization is 
conserved as shown by "1" of Fig. [TJ). The inclination 
breaks the conservation and causes oscillation as shown 
by "2" to "4" of Fig. [Tb. The oscillation would have 
the largest amplitude if 8 — ir/2. Meanwhile, Pmag 2 (t) 
oscillates between zLM 1 /N. It implies that, when the ini- 
tial magnetization is larger, the amplitude is larger. The 



overlap of "4" and "5" in lb implies that the interaction 
plays no role in the oscillation of magnetization. 

(v) A larger uj would reduce the size of the system. 
Since the strength G depends on the spatial wave func- 
tion and would become larger if the size is smaller, the 
effect of interaction would become stronger if u) is larger. 

(vi) Although the evolution caused by the inclined field 
does not depend on N as mentioned, the evolution caused 
by interaction does. Since all the solid curves of Fig. [T^, 
appear to be strictly confined by the dashed curve, how 
these curves vary with A" depends on how the dashed 
curve varies with N . The latter is referred to previous 
literatures [Toj. 

(vii) When the quadratic Zeeman term is taken into ac- 
count, there is no analytical solution. However, the evo- 
lution can be solved numerically We found that, when 
the field is weak (say, B < 0.1 mG), the effect of the 
quadratic term is negligible. 

In summary, an analytical approach beyond the mean 
field theory has been proposed to describe the spin- 
evolution of small condensates. A magnetic field is ap- 
plied along the Z'-axis which is deviated from the Z-axis 
of quantization. Under the single-mode-approximation, 
exact time-dependent solution of the Hamiltonian has 
been obtained. A formula governing the evolution has 
been derived, and related numerical results have been 
presented. It was found that a very weak magnetic field 
can cause a strong, swift, and nearly periodic oscillation 
which is highly sensitive to the magnitude and direction 
of the magnetic field. The high sensitivity implies that 
the phenomenon might be useful for measuring the di- 
rection of a very weak field. When the magnetic field is 
stronger and the inclined angle is larger (say, B is in the 
order of mG and 8 > tt/15), the oscillation cycles will be 
very dense and the observation might be misunderstood 
as random fluctuations. 



Acknowledgments 

This work is supported by the NSFC under the grants 
10874249 and from the project of National Basic Re- 
search Program of China (2007CB935500). 



[1] D.M. Stamper-Kurn, M.R. Andrews, A. P. Chikkatur, S. 
Inouye, H.-J. Miesner, J. Stenger, and W. Ketterle, Phys. 
Rev. Lett. 80, 2027 (1998). 

[2] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H. -J. Mies- 
ner, A. P. Chikkatur, and W. Ketterle, Nature (London) 
396, 345 (1998). 

[3] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[4] T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 
(1998). 

[5] A. Gorlitz, T. L. Gustavson, A. E. Leanhardt, R. Low, 
A. P. Chikkatur, S. Gupta, S. Inouye, D. E. Pritchard, 
and W. Ketterle, Phys. Rev. Lett. 90, 090401 (2003). 



[6] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005). 

[7] A. Sorensen, L.-M. Duan, J.I. Cirac, and P. Zoller, Na- 
ture 409, 63 (2001). 

[8] M.-S. Chang, CD. Hamley, M.D. Barrett, J.A. Sauer, 
K.M. Fortier, W.Zhang, L. You, and M.S. Chapman, 
Phys. Rev. Lett. 92, 140403 (2004) 

[9] M.-S. Chang, Q. Qin, W.X. Zhang, L. You, and M.S. 
Chapman, Nature Physics (London) 1, 111 (2005). 
[10] C.K. Law, H. Pu, and N.P. Bigelow, Phys. Rev. Lett. 81, 
5257 (1998). 

[11] H. Pu, C.K. Law, S. Raghavan, J.H. Eberly, and N.P. 



■5 



Bigelow, Rhys. Rev. A. 60, 1463 (1999). 
[12] L. Chang, Q. Zhai, R. Lu, and L. You, Phys. Rev. Lett. 

99, 080402 (2007). 
[13] S. Uchino, T. Otsuka, and M. Ueda, Phys. Rev. A 78, 

023609 (2008). 

[14] X. Cui, Y. Wang, and F. Zhou, Phys. Rev. A 78, 

050701(R) (2008). 
[15] J. Cheng, H. Jing, and Y. J. Yan, Phys. Rev. A 77, 

061604(R) (2008). 
[16] R. M. Bradley, J. E. Bernard, and L. D. Carr, Phys. Rev. 

A 77, 033622 (2008). 
[17] J. Kronjager, K. Sengstock, K. Bongs, New J. Phys. 10, 

045028 (2008). 

[18] T. Kuwamoto, K. Araki, T. Eno, and T. Hirano, Phys. 

Rev. A 69, 063604 (2004). 
[19] J. Mur-Petit, M. Guilleumas, A. Polls, A. Sanpera, M. 

Lewenstein, K. Bongs and K. Sengstock, Phys. Rev. A 

73, 013629 (2006). 



[20] J. Kronjager, C. Becker, P. Navez, K. Bongs, and K. 

Sengstock, Phys. Rev. Lett. 97, 110404 (2006). 
[21] H. Pu, S. Raghavan, and N.P. Bigelow, Rhys. Rev. A., 

61, 023602 (2000). 
[22] S. Yi and H. Pu, Phys. Rev. A 73, 023602 (2006). 
[23] E. G. M. van Kempen, S. J. J. M. F. Kokkelmans, D. J. 

Heinzen, and B. J. Verhaar, Phys. Rev. Lett. 88, 093201 

(2002). 

[24] W. Zhang, D.L. Zhou, M.-S. Chang, M.S. Chapman, and 

L. You, Phys. Rev. A 72, 013602 (2005). 
[25] Y. Wu, Phys. Rev. A 54, 4534 (1996). 
[26] M. Luo, C.G. Bao and Z.B. Li, J. Phys. B: At. Mol. Opt. 

Phys. 41, 245301 (2008). 
[27] C.G. Bao, Acta Sci. Nat. Univ. Sunyatseni 43, 70 (2004). 
[28] C.G. Bao and Z.B. Li, Phys. Rev. A 72, 043614 (2005). 
[29] M. Luo, C.G. Bao and Z.B. Li, Phys. Rev. A 77, 043625 

(2008). 



